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Magnetic properties of iron chalcogenide superconducting materials are investigated using density 
functional calculations. We find the stability of magnetic phases is very sensitive to the height of 
chalcogen species from the Fe plane: while FeTe with optimized Te height has the double-stripe- 
type (tt, 0) magnetic ordering, the single-stripe-type (tt, tt) ordering becomes the ground state phase 
when Te height is lowered below a critical value by, e.g., Se doping. This behavior is understood 
by opposite Te-height dependences of the superexchange interaction and a longer-range magnetic 
interaction mediated by itinerant electrons. We also demonstrate a linear temperature dependence of 
the macroscopic magnetic susceptibility in the single-stripe phase in contrast to a constant behavior 
in the double-stripe phase. Our findings provide a comprehensive and unified view to understand 
the magnetism in FeSe^Tei-a; and iron pnictide superconductors. 
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In a short period of time since the recent discover- 
ies P, 0, 3 , the research of iron-based superconductors 
(SC) has been extended to a large variety of materials. 
Among them, iron chalcogenides FeSea;Tei_a; are unique 
with their structural simplicity that they are composed of 
only iron-based layers while maintaining the same nomi- 
nal Fe^+ charge state as the iron pnictides. Unlike most 
of iron SCs which have to be substantially doped to sup- 
press the inherent antiferromagnetism (AFM) and de- 
velop the superconductivity, undoped FeSe is not mag- 
netically ordered, superconducting with ~ 9K 0, 
while undoped FeTe is not superconducting but magneti- 
cally ordered. Very uniquely, FeTe has an AFM ordering 
pattern^ so called "double stripe" with (7r,0) ordering 
vector [o, 0], while all other iron-based superconducting 
materials exhibit "single stripe" type ordering pattern 
with the (tt, tt) ordering vector. 

Nature of magnetism in iron SCs has been a controver- 
sial issue between Fermi surface (FS) nesting and local 
spin moment interactions. Several recent works suggest 
unified pictures based on local moment interactions in- 
cluding some role of itinerant electrons but not in the way 
as in FS nesting scenario d, H E [HI, [H, [H, Q El El ■ 
Since FeTe has the (tt, tt) FS nesting similarly to other 
iron SCs ^ITj, the observed (7r,0) magnetic ordering in 
FeTe may indicate irrelevance of FS nesting mechanism 
at least for FeTe [3]. A recent first-principles study [l8| 
suggests the emergence of (tt, 0) FS nesting assuming sub- 
stantial electron doping by excess Fe, however, this dop- 
ing effect was not observed in an angle-resolved photoe- 
mission experiment [l9[. 

The magnetic stability of (tt, 0) ordering over (tt, tt) in 
FeTe can be effectively described by the nearest, second 
nearest, and third nearest neighbor exchange parameters, 
Ji, J2, and J3, respectively, with the condition J3 > J2/2 
However, what makes such a condition hold uniquely 
for FeTe, neither FeSc nor iron pnictides, is not uncov- 



ered. Furthermore, the absence of linear temperature (T) 
dependence of magnetic susceptibility in FeTe @, E3| , in 
contrast to other iron SCs, is yet to be understood. 

In this paper, we perform the first-principles calcula- 
tions to investigate magnetic properties of iron chalco- 
genides. We find that Te height from the Fe plane is 
a key factor which determines AFM ordering patterns in 
FeTe, so that the ground state magnetic ordering changes 
from the (tt, 0) with the optimized Te height to the (tt, tt) 
patterns when Te height is lowered. We observe the 
same effect with FeSe, concluding that it is not chalco- 
gen species but chalcogen position that determines the 
magnetic ordering in FeSea;Tei_j^, and that this feature 
originates from different chalcogen-height dependences of 
Ji, J2, and J3. Our calculated macroscopic magnetic 
susceptibilities (xm) for (tt, tt) and (tt, 0) AFM orderings 
show linear and constant T dependences, respectively, 
suggesting a clue to understand the puzzling T depen- 
dences of Xm measured in iron SCs. Our results, alto- 
gether, provide a comprehensive view on the magnetism 
in FeScxTci-a; and iron pnictide SCs. 

Our first-principles calculations are based on the 
density-functional theory with the generalized gradient 
approximation [iol and the ab-initio norm-conserving 
pseudopotentials as implemented in SIESTA code [21 1. 
Semicore pseudopotentials are used for Fe, and electronic 
wave functions are expanded with localized pseudoatomic 
orbitals (double zeta polarization basis set), with the cut- 
off energy of 500 Ry for real space mesh. Brillouin zone 
integration is performed with 6x6x6 k-point grid us- 
ing a 2a X 2a X c supercell which contains 8 Fe and 8 Te 
atoms, and is commensurate with both the double- and 
single-stripes AFM ordering patterns. 

The atomic structure of FeTe is shown in Fig. 1(a). 
Fe atoms form a plane and Te atoms are bonded to Fe 
atoms above and below the Fe plane, and the distance of 
a Te atom from the Fe plane is defined as the Te height 



2 




FIG. 1: (Color online), (a) Atomic structure of FeTe where 
ZTe is defined with respect to the Fe plane marked by a shaded 
square. Spin densities of (b) AFMl, (c) AFM2, and (d) AFM3 
phases are displayed on the z = plane, and the supercell is 
denoted by a black square in each figure. 



{zTe) as shown in the figure. Structural optimization 
is performed with the unit-cell lattice constants fixed to 
the experimental high-temperature tetragonal structure 
measured at 80 K, a = c = 3.812 A and c = 6.252 A 
0, although the low temperature AFM phase is associ- 
ated with a monoclinic structure. This is for simplicity 
in defining ZTe and easy comparison with FeSe case. We 
have checked that this simplification does not alter our 
main conclusion. We consider ferromagnetic (FM) phase 
and three different AFM phases to study the magnetic 
properties: 'checkerboard' (Fig. 1(b)), 'single stripe' 
(Fig. 1(c)), and 'double stripe' (Fig. 1(d)) type order- 
ings, which hereafter are referred to as AFMl, AFM2, 
and AFM3, respectively. 

To investigate effects of chalcogen height on magnetic 
orderings, we calculate energetic stabilities of magnetic 
phases as a function of zte- Starting from the optimized 
ZTe of 1.81 A, the total energies of FM, AFMl, AFM2, 
and AFM3 phases are calculated with decreasing ztc as 
shown in Fig. 2(a). At optimized ztci AFM3 is most 
stable in agreement with the experiment and the pre- 
vious calculation 0] , and AFM2 is slightly higher in 
energy by 32 meV per Fe. Interestingly, FM phase is 
found to be more stable than AFMl phase in contrast to 
our previous experience on iron pnictides that FM phase 
relaxes to the non-magnetic (NM) state without a con- 
straint on the magnetic moment and has higher energies 
than AFM phases even with a fixed spin moment (FSM) 
calculation. When ztc is decreased to 1.75 A, AFM3 is 
still more stable than AFM2 but the energy difference 
reduces. After ZTe is further decreased approximately 
below 1.71 A (see the vertical arrow in Fig. 2(a)), finally 
AFM2 phase turns more stable than AFM3 phase, and 
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FIG. 2: (Color online), (a) Total energy versus zte calculated 
for FM, AFMl, AFM2, and AFM3 phases, (b) Exchange pa- 
rameters J's are shown as functions of zxe which are eval- 
uated with fixed Fe magnetic moment of 2.95 /is. Vertical 
arrows denote the critical ZTe's where AFM2-AFM3 switch 
takes place which are slightly different in (a) and (b) due to 
different treatment of local magnetic moments (see text). 



the relative stability of AFM2 over AFM3 phases is even 
enlarged as ztc is decreased further to 1.63 A, towards 
our optimized Se height 1.55 A of FeSe. The Fe magnetic 
moment is decreasing along with ZTe, and the averaged 
values over the three AFM phases for each zxe are 2.94, 
2.84, 2.72, and 2.58 ^lB for zxe = 1-81, 1.75, 1.69, and 
1.63 A, respectively. Except for the ztc =1-81 A case, we 
need to perform the FSM calculation to obtain the FM 
phase, in which the Fe moment is fixed to the average 
value of the three AFM phases with the same zxe- 

Since AFM3 ordering in FeTe is apparently incompati- 
ble with the FS nesting scenario, we consider Heisenberg 
type local moment interactions, and estimate the three 
exchange parameters Ji , J2, and J3 , as functions of ztc , 
from the calculated total energies of each magnetic phase 
and ZTe- Since J's depend on the magnitude of the mag- 
netic moment in general (2^ . we perform the FSM cal- 
culations, with the Fe spin moment being fixed to 2.95 
/is (the value for AFM3 at ZTe=l-81 A), for all the mag- 
netic phases and ztc- The result is shown in Fig. 2(b). 
At optimized ztc — 1-81 A, J2 is largest, while Ji is es- 
timated to be negative as the FM phase is more stable 
than AFMl phase at this ZTe, and J3 is almost compa- 
rable with J2. As ZTe decreases, we see a clear tendency 
that J I and J2 increase while J3 decreases. Around ztc 
= 1.72 A, J2/2 starts to exceed J3 which reflects that 
AFM2 phase becomes more stable than AFM3 from this 
point on. We should note that this critical value of ZTe 
is slightly different from the one in Fig. 2(a) where the 
magnetic moments are unrestricted, while they are fixed 
in the estimation of J's as mentioned above. 

In a word, the switch of ground state magnetic phase 
between AFM2 and AFM3 with varying ztc is due to 
opposite ZTe dependence between J2 (along with Ji) and 
J3. The behavior of Ji and J2 can be understood rela- 
tively easily if we recall their superexchange nature [23 | . 
As Te atom gets close to the Fe plane, Fe-Te-Fe angle 
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FIG. 3: (Color online). Magnetic susceptibility XaiXyr' = 
Q\ui = 0) of the NM phase is shown for (a) zte^I.SI and 

(b) ZTe = 1.63 A. Delta- like magnetic field is applied at the 
center of the figure. PDOS in the NM phase is displayed for 

(c) Fe 3d and (d) Te 5p states. 



approaches 180° which maximizes the overlap between 
Te 5p and Fe Zd orbital robes, and the superexchange in- 
teraction should get stronger. As for J3, a longer-range 
interaction is needed where isolated local spin moments 
can be coupled via itinerant electrons around the Fermi 
energy {Ep), such as so-called RKKY interaction. 

To understand ztc dependence of J3 , we consider how a 
localized magnetic perturbation polarizes the surround- 
ing initially spin-unpolarized electron cloud. We apply 
a delta- like magnetic field on a Fe atom of the FeTe 
system in the NM phase, by adding a potential differ- 
ence between spin up and down electrons at the center of 
the Fe atom, and obtaining self-consistent electron spin 
density in response to this external field, as displayed in 
Figs. 3(a) and (b). By definition, this plot corresponds to 
Xs(r, r' = 0; a; = 0) when the delta field is applied to the 
origin. Here we use 2x2x1 enlarged supercell of our orig- 
inal supercell shown in Fig. 1 to reduce the effect of other 
external fields from the periodic images of the supercell. 
When Te is at the optimized height 1.81 A (Fig. 3(a)), 
the electron cloud is strongly spin-polarized by the delta 
field and show an alternating pattern, with a slow decay 
rate with distance from the origin. In the meanwhile, 
when ZTe is lowered down to 1.63 A in Fig. 3(b), the 
spin polarization decays much more rapidly away from 
the origin. Hence, Xs(r, r' = 0;w = 0) for the NM phase 
becomes much more short-ranged with smaller zxe , con- 
sistently with the fact J3 decreases with decreasing ztc- 

The ZTe dependence of Xs can be understood by ana- 
lyzing the change in the electronic structure along with 



FIG. 4: (Color online), (a) Macroscopic magnetic suscepti- 
bility xm versus electron temperature is shown for two AFM 
orderings. (b) DOS of FeTe for AFM2 and AFM3 phases. 



varying zxe- Figure 3(c) represents the projected density 
of states (PDOS) on Fe 3d, where two peaks around -3 
eV and Ep have t2g character. When Te is lowered, the 
t2g states couple to Te 5p states more strongly so that the 
corresponding peaks around -3 eV in Fig. 3(c) and (d) 
move down to the higher binding energy, gaining more 
Fe t2g weight and losing Te 5p weight with an effective 
charge transfer from Te to Fe. Then the density of states 
(DOS) at Ep is greatly reduced due to the Fe t2g weight 
transfer to the states around -3 eV and overall charge 
transfer to Fe states. This in turn results in generally re- 
duced magnetic susceptibility Xs, since it depends on the 
number of states around Ep (which correspond to 'itiner- 
ant electrons'). The higher DOS at Ep for high zxe also 
explains the stabilization of the FM phase at ZTe=1.81 
A shown in Fig. 2, according to the Stoner criterion. 

The Curie- Weiss-like behavior, i.e., the 1/T depen- 
dence of magnetic susceptibility is a peculiar feature mea- 
sured in FeTe 0, [l^l which undoubtedly shows the lo- 
cal nature of magnetism, while other iron-base SCs show 
linear-T dependence. Motivated by the fact the ground- 
state magnetic phase is AFM3 for FeTe and AFM2 for 
other iron SCs, we evaluate the macroscopic magnetic 
susceptibility xm for different magnetic orderings, as a 
function of electron temperature by applying a uniform 
magnetic field and obtaining self-consistent electron spin 
density. The result is shown in Fig. 4(a). The difference 
is obvious: xm depends linearly on T for AFM2 while it 
is almost constant for the AFM3 phase. This difference 
results from the dramatic difference in DOS at Ep, shown 
in Fig. 4(b). AFM2 phase has a large dip in DOS at Ep 
because of the band-repulsion between electron and hole 
bands which are separated by the q = (tt, tt) vector in the 
Brillouin zone of the NM phase [ll], while DOS at Ep 
is large in AFM3 phase because the perturbing potential 
by the magnetic ordering has only q = (tt, 0) component 
which cannot couple the electron and hole FSs. Since 
the susceptibility is larger when there is large number 
of states around Ep, the overall value of xm is smaller 
for AFM2 phase. However, as the electron temperature 
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increases with the Fermi-Dirac distribution, more states 
become available to participate in the magnetic response 
and xm rises. 

For Xm in high-temperature paramagnetic phase, we 
consider canonical ensemble of local spin moments, in 
which the probability of each configuration is given by 
the Boltzmann factor associated mainly with the Heisen- 
berg Hamiltonian. When a uniform magnetic field is ap- 
plied, Xm has 1/T-like dependence if the spin moments 
are frozen in each configuration. With itinerant electrons, 
Xm also has contribution from field-induced spin-density 
change in each configuration averaged by the Boltzmann 
factor, which is dominated by the ground state and low- 
energy configurations possibly composed of disordered 
domains of the ground-state AFM phase [2^ . Then, the 
spin-density change in the ground-state AFM will con- 
tribute significantly to xm in the paramagnetic phase. 
Thus we suggest xm in the paramagnetic phase has a 
linear-T component for lower chalcogen height, while the 
component is absent for higher chalcogen height. 

We also have checked that the ground states in FeSe 
and LaFeAsO, which are AFM2 with optimized atomic 
positions, change to AFM3 when Se or As height is in- 
creased. This implies the magnetism in FeTe in fact is 
not qualitatively different from that of other iron-based 
SCs, and the peculiar (tt, 0) ordering pattern is just a spe- 
cific realization among many possible magnetic orderings 
determined by the relative strength of local-moment ex- 
change interactions J's, which can be controlled contin- 
uously by varying the chalcogen (and pnictogen) height, 
for example, by alloying such as FeSe^jTei-x. Although 
FeSe and other alloys with FeTe are not experimentally 
reported to exhibit any magnetically ordered phase pos- 
sibly due to the off-stoichiometry etc., our results sug- 
gest that AFM2 is the lowest-energy configuration in the 
canonical ensemble for paramagnetic FeSej;Tci_2, with 
X over a critical value. This possibly results in the 
linear-T dependence oixM as demonstrated above. How- 
ever, when Se concentration is not enough, the chalcogen 
height is so high that AFM3 is the lowest-energy con- 
figuration, and the Curie- Weiss-like 1/T dependence of 
Xm @, [3 , ordinary for the local moment magnetism, is 
observed in the paramagnetic phase. 

Increasing Se in FeSea;Tei_a; is reported to suppress 
(tt -f 6, 0) spin fluctuation and enhance (tt, tt) fluctuation 
with the superconductivity ^2^. Our result shows that 
this might be related to the continuous change of the rela- 
tive energy between AFM3 and AFM2 orderings with in- 
creasing X. When x is small and chalcogen height is large, 
the AFM3 state and (tt, 0) spin fluctuation has greater 
probability in the canonical ensemble for the paramag- 
netic phase, weakening the interaction between electron 
and hole FSs separated by g = (tt, tt). As a; increases over 
the critical value, the AFM2 state and (tt, tt) spin fluctua- 
tion start to take over, which coincides with the supercon- 
ductivity observed in FeSe^Tei-^:. Therefore, our result 



shows that the superconductivity in iron chalcogenides 
may be related to (tt, tt) spin ordering and fluctuations 
as likely as in iron pnictides. Our findings further sug- 
gest that the superconducting state might be realized in 
FeTe by controlling the Tc height by, such as, applying a 
biaxial strain. 

In conclusion, we show that the unique magnetic or- 
dering in FeTe is due to the relatively high chalcogen 
height compared with FeSe as the underlying magnetic 
interactions depend on the chalcogen height in different 
ways. This behavior is found to apply for other iron 
chalcogenide and iron pnictides, hence the magnetism in 
iron-based SCs including FeTe can be understood with 
a same unified mechanism. Temperature dependence 
of magnetic susceptibility evaluated for the AFM2 and 
AFM3 phases sheds light on the puzzling experimental 
results observed for FeSe^jTei-a:. 
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